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Abstract A new method of virtual unknown parameter is proposed to synchronize 
two different systems with unknown parameters and disturbance in finite time. Virtual 
unknown parameters are introduced in order to avoid the unknown parameters from 
appearing in the controllers and parameters update laws when the adaptive control 
method is applied. A single virtual unknown parameter is used in the design of 
adaptive controllers and parameters update laws if the Lipschitz constant on the 
nonlinear function can be found, while multiple virtual unknown parameters are 
adopted if the Lipschitz constant cannot be determined. Numerical simulations show 
that the present method does make the two different chaotic systems synchronize in 
finite time. 
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1. Introduction 

Synchronization in chaotic dynamic systems has attracted increasing attention of 
scientists from various research fields for its advantages in practical applications [1]. 
A wide variety of methods have been proposed for synchronization of chaotic systems, 
including linear feedback control [2,3], sliding mode control [4,5], adaptive control 
[6,7] and so on. Most of the methods mentioned above are used to guarantee the 
asymptotic stability of chaotic systems. In other words, the convergence of the 
synchronization procedure is asymptotic with infinite settling time. 

In many practical applications, one of the main concerns is that the behavior of 
system is stable over a finite-time interval. To achieve faster convergence in control 
systems, finite-time control is a very useful technique. The finite-time control 
techniques have demonstrated better robustness and disturbance rejection properties 
[8-10]. For the nonlinear system of the form [11] 

l*j =x"*l l +0M,---x J ) (7 = 1, - "-I) 
[x n =u + <t> n ( Xl ,---x n ) 

where rrij are odd integers, output feedback controllers can be explicitly constructed 

to globally stabilize the system in finite time [11-15]. Finite-time observers for 
nonlinear systems that can be put into a linear canonical form up to the output 
injection have been proposed [16-18]. Whereas, linearization method just ensure the 
local stability. Methods in [11-15] are suitable only for the special systems mentioned 
above, and cannot be applied directly to synchronization for general chaotic systems, 
such as Lorenz system and Chen system. 

So far, adopting the concept of convergence in finite time [19], finite-time 
sliding-mode controllers and feedback controllers have become two effective methods 
in finite-time synchronization. The state feedback control law involving the term 

-k\Sign{x) |xp - k 2 sign(v)\v\ a2 is developed to solve the finite-time synchronization 

problem of two different second-order chaotic systems without uncertainty [8]. 
However the controller in [8] contains all the information appeared in the error 
dynamical system or the drive-response systems, which is hard to be realized in 
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practical application. 

Nevertheless uncertainties do exist in any real world systems. The finite-time 
control problems for linear systems subject to time -varying parametric uncertainties 
and to exogenous constant disturbances is considered[20]. A sliding-mode controller 
is designed to drive a class of chaotic systems with nonlinear inputs to track a smooth 
target trajectory in finite time[21]. The nonsingular terminal sliding mode control is 
used to stabilize the chaotic systems with uncertain parameters or disturbances in 

finite time [22]. A controller containing the term e p with /? = q I p < 1 is proposed to 

realize the finite-time synchronization of two identical unified chaotic systems with 
uncertain parameters[23]. In view of a control Lyapunov function, controllers 
including the term of Lie derivative (L B V) T L f V are designed to cope with the chaos 

synchronization for the coupled identical united chaotic system with uncertain 
parameters [24]. 

In this paper, we introduce a virtual unknown parameter method to realize the 
finite-time synchronization between two different systems with uncertainties. In order 
to avoid the unknown parameters from appearing in controllers and parameters update 
laws, virtual unknown parameters are introduced. With the help of a virtual parameter, 
an adaptive controller and the corresponding parameters update laws are therefore 
designed to synchronize two coupled chaotic systems with uncertainties in finite time 
when the Lipschitz constant on the nonlinear function can be found. Sometimes it is 
hard to determine the Lipschitz constant. In order to overcome this difficulty, we 
introduce multiple virtual unknown parameters to develop techniques for obtaining 
proper adaptive controllers. Some typical chaotic systems, Rossler system, Chua's 
circuit and Lorenz system, are chosen as illustrative examples to demonstrate the 
effectiveness of the proposed method. 

2. Problem statement 

Consider a class of chaotic systems with unknown parameters and disturbance in 
the form of 
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(i) 



where xeR n is the state vector , a e R m is an unknown parameter vector, 
/(x)ei?° xl and F(x) e i?" xm are nonlinear matrix functions, d{x,t)&R" is the 

external disturbance. Let Q t c i?" be a bounded region containing the whole 

attractor of drive system (1) such that no trajectory of system (1) ever leaves it. This 
assumption is simply based on the bounded property of chaotic attractor. Also, let 
M a R m be the set of parameter under which system (1) is in a chaotic state. 
Introducing a virtual unknown parameter a m+l , Eq.(l) can be transformed into 



x = f(x) + F(x)a + d(x,t) 



(2) 



where /(x)ei?" xl and F(x)eR 



nx(m+l) 



are nonlinear matrix functions , 



a T = (a T ,a m+l ) e R n+l is a new parameter vector, and the term in F(x) 
corresponding to virtual unknown parameter a m+l is called virtual term. 

Many chaotic systems can be transformed into the form of system (1) and (2). Let 
us take Lorenz system as an example, 



A, 



\ X 3j 



^ (x 2 -x x OYtf^ 
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V X l X 2 J 







x x 



Ox 



3 7 



= f(x) + F(x)a 



V X l X 2 J 



x 2 — Xj 







x -x 



x 3 



a-, 



a, 



= f(x) + F(x)a. 



where unknown parameters, a 4 = 1 is a virtual unknown parameter, 

and -x 2 is the virtual term corresponding to virtual unknown parameter a 4 . 
Remark 1. Each component, linear or nonlinear, of f(x) can be elected as a virtual 



term in F(x) corresponding to parameter a m+l , where a m+l ^ . 



3. Finite-time chaos synchronization 

3.1 A single virtual unknown parameter in response system 

Let Eq. (2) be the drive system. The corresponding response system is another 
different system defined as below, 

y= s( i+ a M + < % (3) 

which can be transformed into 

y = g(y) + G(yW + u(t,x) (4) 

where y^R" is the state vector, g(y) e R"* , g(y) e R nxl ,G(y)^R nx ' and 

G(y)eR nx{M) are nonlinear matrix functions. W = (P T ,f3 M ) e R M is a new 

parameter vector, in which /? e i?' is an unknown parameter vector of the initial 

system and J3 M is a virtual unknown parameter. u(t,x) is a controller. Let 

Q v ci?" be a bounded region containing the whole attractor of response system (4) 

with u(t,x) = 0. 

The object of this paper is that the coupled system (2) and (4) can be synchronized 
in finite time by designing an effective controller u(t,x) and adaptive laws. In order 
to do so, some hypothesizes are made as follows, 

HI. The nonlinear functions /(•) and g(-) are smooth and continuous on a 
bounded closed region Q containing both Q x and Q y . Thus there exists a 
constant M such that 

\f{xy- g {\% i. 

H2. The nonlinear matrix function g(-) satisfies Lipschitz condition, namely 

||£(*)-,gO0||^||*-;HI 

where L- is an appropriate positive constant. In this paper, ||*|| denotes matrix or 

m n 

vector norms, defined as ||^|| = ^ a\ )' 12 for the matrix A = (a !7 )„ x ,„ . 
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Fortunately, many chaotic systems transformed into the form of Eq.(2) and (4) 
satisfy assumptions HI and H2, such as Lorenz system, Rossler system, Genesio-Tesi 
system. 

H3. The uncertain parameters a, ft and disturbance d(x, t) are all norm bounded, 

\\ a \\<s a \p\\<s p ,\\d(x,t)\<s d 

where 8 a ,8p,8 d are positive constants. Then there exist two positive constants S a 

and 5 p , such that 

||«|| = 5 a <8 a + \a m+x | , = S p <S p + \fi M | . 
Introduce the following notations, 

which are m + 1 -dimensional and / + 1 -dimensional constant vectors , respectively. 

Define the error variable as e—x—y. Subtracting (4) from (2) yields the 

synchronization error dynamics as 

e= A $--f£ $~ Eft*' drx, 1) (5) 

So the problem to realize synchronization between systems (2) and (4) in finite 
time is transformed into the finite time stability of error dynamics (5). Hence the 
definition of finite-time synchronization and some necessary lemmas are introduced 
in the follows. 

Definition [23] Consider the following two chaotic systems: 

X m = f( X m )' ^ 

x s =h(x m ,x s ), 

where x m ,x s are two n -dimensional vectors. The subscripts ' m ' and ' n ' stand for 

the master and slave systems, respectively. / :R n — > R" and h:R n ^> R" are 
vector-valued functions. If there exists a constant T > , such that 

and |x m -xj = 0, if t>T , then synchronization of the system (6) is achieved in a 
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finite-time. 

Lemma 1[23] Assume that a continuous, positive-definite function V(t) satisfies 
the following differential inequality, 

V{t)< - cV(, W>0,F(O>0, 
where c>0,0<^<1 are all constants. Then, for any given t , V(t) satisfies the 
following inequality, 

(t) < v 1 -" (t ) - c(l - tjXt -t ),t <t< tl , 



and V(t) = for t>t l with t 1 =t + 



cG-17) 

Lemma 2[1 1] For any given a,b e R, < g < 1 , the following inequality holds, 

(\a\ + \b\) q <\a\ q +\b\ q . 
Now we can establish the following result. 

Theorem 1. The coupled system (2) and (4) can be synchronized in finite time if 
the hypotheses H1-H3 hold and the following conditions (I) -(IT) are satisfied, 

(I) The controller is chosen as 



u(x, t) = k x e + (k 2 + 1) TT-7 + F(x)a - G(y) f3 



IK 
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(7) 



+2 



where k x > L g ,k 2 > S d + M , a and /? are two vector parameters. 

(n) The update laws of the parameters a and (i are 

A„ -a 



a =F(x) r e + - 



P = -G(yfe + 



(8) 
(9) 



where d m+l and represent the (m + Y)-th and (/ + l)-?/z component of 



vectors a and /J , respectively. 

For a fixed t , the finite time t { is determined by 

t x =t Q + 2W^(t ) . (10) 
where V is a positive-definite function satisfying Lemma 1 . 
Proof: Choose a Lyapunov function of the form 

Consequently, the time-derivative of V along the error system (5) is 
V = e T e-(a- a) T a - (fi - P) T p 
= e T (fix) - g(y) + F(x)a - G(y) j3 + d(x, t) -u(t, x)) -(a-af h-(P - Pf P 
From above hypothesizes H1-H3 and Eq.(7), one can obtain 

V = e (f(x) - g(x)) + e T (g(x) - g(y)) + e T F(x)(d -a)- e T G(y)0 - P) 



+e T d(x,t)-e r [k l e + (k 2 +l) — + 2 



-ia-af&-iP-P) T P 



+ A„ • \\d\\ 



- + 2- 





|2 






P 




P M - 


~P M 





< M \e\ + L \\e\\ + e T F(x)(a -a)- e T G(y)(j3 ~P) + 8 d \e\ 



\K 


2 

+ 




•||« 











■ + 2- 



-e T [k l e + (k 2 + l)-^- + 2^ 



-ia-af&-iP-P) T P 
By using k x >Lk 2 >5 d +M and Eqs.( 8) and (9), it yields that 





2 1 
+ 






P 




p M - 


~P M 





V<-\\e\\-e ^M + \K\\fl 



- + 2- 



-(a-d) 



A — oc — " t A„-jS 

' ■ -(P-py " ' 



Using the inequalities 



\a m+l -a m+l \<\\a-a\\ and 



p M -p M < p-p 





|2 

















r] 



(11) 



we have 
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1 



< - 



1 



and 



\ a m+\ a m+\\ 



\\a -a\\ 



1 

P M -P M 



< 



1 

p-p 



(12) 



With the help of inequality (12), inequality (11) can be transfonned into 

+ (a-d) 



''< H (« a)'-! 1 a)'- ° ! \ (p /i)' p " ' P 



\\a — a\\ 



«» + l-«m + l| 



P-P 



HP-P) 

P M ~P M 
The following expressions hold, 



|A || 2 +||A 

I a || || a 


• ||a|| e 






ii ii 2 



+ 2- 
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+ 






p 






-A + i 





(13) 



(a - a/ • (a - A B ) = llalf - a 7 A a - a r a + a r A a < 2 ||A a If + 2 ||A 



\\a\\ 



(P-py-(P-& p )= 

This leads to 

F<-||e||-||«-«||- 
From Lemma 2, we obtain 



■P T A S -P T P + P T A B <2\A B \ +2 AJ- ^ 



_ .2 



K<-(e + a -a + ) 2 



i i 

= _ 2 2J72_ 



By applying Lemma 1, it can be concluded that the error system (5) is stabilized at 
origin point in finite time t x defined by 



t,=t n + 



V 2 (t ) 



o 1 1 



= < + 2=K ! (O 



2>0-i) 

That is, the response system (4) synchronizes the drive system (2) in finite time. 
Remark 2. We try to construct a suitable controller such that the specific Lyapunov 
function V can satisfy Lemma 1 which guarantees the finite-time synchronization. 

The unknown parameter vectors a and P in V should not appear in controller and 

parameter update laws in practice. Using the method, for example, in [25], the 
controllers contain the unknown parameters. If we use the same method as above but 
do not use the virtual unknown parameters, the controller and parameter update laws, 
i.e. Eq.(7)-(9), should be 

9 



u(x,t) = k l e + (k 2 + Y) — + F(x)d-G(y)/3 



A + A -a e 
+2 \\ aJL II « II H ._g^ + 2 j 

Fi-«» llell 



^ = F(x) r e+ A " a . , /? = -G(j/) r e + 
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2 

+ 






ft 


e 






|e| 2 



A -A 



which contain unknown parameters, namely a t and /? y . , of the drive and response 

systems, respectively. Theorem 1 indicates that introduction of a virtual unknown 
parameter is available to avoid the unknown parameters from appearing in controllers 
and parameters update laws and realize the finite -time synchronization of two 
different chaotic systems. 

Remark 3. Obviously, the magnitude of — =- in the controller u(x,t) will increase 



6 6 

unboundly and become infinity when e — > . So, we take — = instead of — : 



r +e \\ e \\ 
in practice, where s is a sufficient small positive constant. This idea also can be 
found in [26]. 

3.2 Multiple virtual unknown parameters in response system 



In this part, we consider the finite time synchronization of two different chaotic 

systems with multiple virtual unknown parameters in response system. Sometimes it 

is hard to determine the Lipschitz constant. For example, if we take the Lorenz system 

presented in Section 2 as a response system, we have 

n 

= >/(-*i* 3 + y^f +{x x x 2 -y<y 2 f 



f(x)-f(y) = 



v x \ x 2 ~ y^i 



Apparently, it is difficult to find a constant L satisfying ||/(x) - f(y)^ <L\\x-y\\ . 
Multiple virtual unknown parameters are introduced to overcome this difficulty. The 
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response system with controller u{t,x) is chosen as 

y = h(y) + H(y)y + u(t,x), (14) 
where y e R" is the state vector, h (y) e R" does only contain linear terms and 
H(y)<ER" x(l+s) is a nonlinear matrix function. y T = (y T ,y M --- y l+s ) e R l+S (s <n) 
is a new parameter vector, in which y e R' is an unknown parameter vector of the 
initial system and y M ---y !+s are virtual unknown parameters. u(t,x) is a controller. 
The hypotheses HI and H2 in Theorem 1 are modified as, 

H T . The nonlinear functions / (•) and h (•) are smooth and continuous on a 
bounded and closed region. Thus there exists a constant M, such that 

\\f( X y~h(^ /• 

H 2' . The nonlinear matrix function h (•) satisfies Lipschitz condition, namely 
\\h(x)-h(y)\\<L j; \\x-y\\, 

where is an appropriate positive constant. Then the following result is attained. 

Theorem 2. The coupled system (2) and (14) can be synchronized in finite time if 
the hypotheses HI', H2'and H3 hold and the following conditions (I)-(n) are 
satisfied, 

(I) The controller is chosen as 

u(x, t) = k x e + (k 2 + 1) tt^t + F(x)a - H(y)y 
IN 

„ ||2 || A ||2 || A || ,,,,, (15) 

A + A • \\d\\ p A + A • y » 

+2 ¥M n^A FII ._L. +2 II o II dLEI ._L.. (i</< 5 ) 
|H| \y l+i -y l+i \ \e\ 

where k l >L- R ,k 2 >8 d +M l , a and ^ are two vector parameters, and 
A y = {5 Y ,5 r ,---5 y ) T is a (/ + ^) -dimensional constant vector, where ||^j| = <? 
^S y +\y m \ + ---+\y l+s \,\<s<n. 

n 



(n) The update laws of the parameters a and y are 

&=F{x) T e + - a ~ a - , (16) 

f = -#(z) r g + (l<i<s) (17) 

where d m+1 and ^ /+( . represent the (m + Y) — th and (l + i)—th component of 
vectors a and ^ , respectively. 

i i 

Furthermore, for any given f , the finite time f, is determined by ^ = f + 2 2 F 2 (^ ) , 
where V is a positive-definite function satisfying Lemma 1 . 

So we can draw a conclusion that the response system (14) can globally 
synchronize the drive system (2) in finite time. The proof procedure of Theorem 2 is 
similar to that of Theorem 1 and omitted here. 



4. Numerical simulation 

To verify the feasibility of the proposed method, the Rossler system and modified 
Chua system, Rossler system and Lorenz system are used as examples in the 
simulations. 



Example 1 The Rossler system and modified Chua system for the case of Theorem 1 

The drive system is Rossler system, with disturbance d{x,t) and unknown 
parameter a described by 
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K d 3 (x,t) 



(18) 



The controlled modified Chua' circuit is used for the response system[27] 
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(19) 
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u 2 (x, t) 






V W 3 C^sO y 



In simulation, values of uncertain parameters of the drive and response systems 



are chosen as a x = -\,a 2 = 0.2,a 3 = -l,/3 l = 10, /? 2 



LOO 



The virtual unknown 



parameters are a 4 = fi 3 = 1 . The small positive constant is set to e = 10 -4 , and the 

disturbance is d = (d l ,d 2 ,d 3 ) T = (0.05x 2 smt,0.\x 3 sm3t,0.0lx l cos2t) T . Figs. 1-2 
shows that the drive system (18) is chaotic under the above parameter values and the 
initial value (x 1 (0),x 2 (0),x 3 (0)) r = (1,1, l) r . From Figs. 1-2, the bounds of the 

chaotic attractor are -10 < Xj < 11.5,-11 < x 2 < 8,0 < x 3 < 24 . Therefore, the 
following auxiliary values can be obtained 

||<5U|<2.7, ||£j<7.lJ<?J<8.1, |<y<17.6,|JJ<18.6 



f(x)-g(x) = 







+ 0.2y 



= > /(-jc 1 +jc 2 ) 2 + (jc 1 jc 3 +0.2) 2 < 278, 



It is easy to verify that g(-) satisfies Lipschitz condition, 

( ^ 



If ( *y ~g( t> 



1 * i~x f~y 




||2- |c. 



By Theorem 1, we chose k { = 3 and k 2 = 281 . The controller is given by the following 
expression, 
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u = 3e + 



281 



e + F(x)a-G(y)fi 



+2 KII +IK 



\\oc\\ 
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ii M 2 
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\l-d A 



where a,/? are updated according to the following adaptation algorithm, 



a = 



x 2 +x 3 ^ 
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-3^2 +-(2^-^) 








o y 2 



e + 



^18.6-/?^ 

18.6-A 
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The initial conditions for the response system and update laws of the parameters 
are (^(0), y 2 (0), y 3 (0)) T = (-1,2,-1/ , (a,(0),a 2 (0),^(0),a 4 (0)) r = (1,2,3, 4) T , 

(A( ).A(°).A(°)) r =(0.5,l,1.5) r , respectively. From Fig.3, one can see the error 
vector e converges to zero aW = 0.002 with control operation. This implies that 
the trajectories of the response system converge to those of the drive system in finite 
time. While by the Eq.(10), the converge time is t = 16.8148 . The difference between 
the numerical and theoretical converge time mostly comes from the conservative 

chose of k x and k 2 . 



Example 2 The Rossler system and Lorenz system for the case of Theorem 2 

Let Eq.(17) be a drive system. The controlled Lorenz system with following 
equation is used for response system, 
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75 J 



= AO0+#O0r- (20) 



Apparently, finding a Lipchitz constant for nonlinear function h(y) needs the 

bounds of the response system, which are difficult to be determined in practice. 

Introduction of multiple virtual unknown parameters can overcome this difficulty 

Values of uncertain parameters of the response system are chosen as 
g 

y l = 1 0, y 2 = — , y 3 = 28 , and the virtual unknown parameters are y 4 =y 5 =l. Therefore, 
the following auxiliary values can be got 

||^ I < 30, \8 y \ < 32, A r = (32, 32, 32, 32, 32) r . 



/(*)-*(*) = 





x 2 

v XjX 3 + 0.2 



= *Jx 2 2 +(x l x 3 +0.2) 2 < 276.5, 



It is easy to verify that h (•) satisfies Lipschitz condition, 

' v 
-x 2 +y 2 




\h(x)-h(y)\ = 



-\\ x ~y\\ 



Proceeding as before, we chose k x = 2 and£ 2 = 280 . The controller is given by the 
following expression, 

280 — — 
u = 2e + -7—r e + F(x)d - H(y)y 
M 



+2- 1 



A I + A„ • \\a\\ 



\l-a A 



P A + A • y 



i-r 5 



where y is updated according to the following adaptation algorithm, 
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The initial condition of y is (n(0),7 2 (0),;r3(0),;f 4 (0),;f 5 (0)) r = (5,6,7,0,l) r . The 

other values of parameters and the adaptation laws of a are the same as those of 
Example 1. From Fig. 4, one can see the error vector e converges to zero at 
t = 0.0008 with control operation. By the Eq.(10), the theoretical converge time is 
t = 29.2494. 



5. Conclusion 

This work presents a general method for synchronizing two different chaotic 
systems with disturbance and unknown parameters in finite time. We introduce virtual 
unknown parameters to avoid the unknown parameters from appearing in controllers 
and parameters update laws. A single virtual unknown parameter is used to design an 
adaptive controller and the parameter update laws which guarantee the finite-time 
synchronization between two coupled chaotic systems. For the case that Lipschitz 
constant is hard to be determined, multiple virtual parameters are taken to develop 
techniques for obtaining proper adaptive controllers. Numerical simulations on the 
basis of the Rossler system and Chua's circuit, Rossler system and Lorenz system are 
presented to show the effectiveness of the proposed method. 
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Figure captions 

Fig.l x 1 -x 2 plane projection of Rossler attractor 
Fig.2 x 1 -x 3 plane projection of Rossler attractor 

Fig. 3 Synchronization error between systems (18) and (19) with initial conditions 

(x x ,x 2 ,x 2 ) = (1,1,1) and iy x ,y 2 ,y 3 ) = {-\,2,-\). 

Fig.4 Synchronization error between systems (18) and (20) with initial conditions 

(x l5 x 2 ,x 3 ) = (1,1,1) and iy v y 2 ,y 3 ) = (-1,2,-1). 



19 



5 

=t] (to 



Fig.l x,-x 2 plane projection of Rossler attractor 




Fig.2 Xj - x 3 plane projection of Rossler attractor 
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Fig. 3 Synchronization error between systems (18) and (19) with initial conditions 
(x x ,x 2 ,x 3 ) = (1,1,1) and (y v y 2 ,y 3 ) = (-1,2,-1). 
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Fig. 4 Synchronization error between systems (18) and (20) with initial conditions 
(x 1? x 2 ,x 3 ) = (1,1,1) and (y 1 ,j 2 ,j 3 ) = (-1,2,-1). 
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